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Abstract 

Recently it was shown how to formulate the finite-element equations of mo- 
tion of a non-Abelian gauge theory, by gauging the free lattice difference equa- 
tions, and simultaneously determining the form of the gauge transformations. 
In particular, the gauge-covariant field strength was explicitly constructed, lo- 
cally, in terms of a path ordered product of exponentials (link operators). On 
the other hand, the Dirac and Yang-Mills equations were nonlocal, involving 
sums over the entire prior lattice. Earlier, Matsuyama had proposed a lo- 
cal Dirac equation constructed from just the above-mentioned link operators. 
Here, we show how his scheme, which is closely related to our earlier one, can 
be implemented for a non-Abelian gauge theory. Although both Dirac and 
Yang-Mills equations are now local, the field strength is not. The technique 
is illustrated with a direct calculation of the current anomalies in two and 
four space-time dimensions. Unfortunately, unlike the original finite-element 
proposal, this scheme is in general nonunitary. 



*E-mail: kmilton@uoknor.edu 
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I. INTRODUCTION 



An alternative approach to lattice field theories, based on the finite-element equations of 
motion, has been under development for over a decade. (For a recent review see [jl|.) Shortly 
after the introduction of this method, it was seen how a Abelian gauge field could be coupled 
to a fermion in this way 0. The resulting Dirac equation was nonlocal: In Minkowski space- 
time, the term proportional to 7-' involved a sum over all values of the corresponding lattice 
coordinate, mj, 1 < mj < M, where M is the number of lattice sites in the j direction, 
while the term proportional to 7° involved a sum over all previous times, < n' < n, where 
n is the current lattice time. Shortly after our paper appeared, Matsuyama |§| proposed a 
local finite-element Dirac equation for QED, based on the immediate introduction of link 
operators into the free Dirac equation. Although it could be argued that this latter approach 
was somewhat unnatural because it introduced interactions into the mass term, the primary 
reason this idea was not pursued was that it was quite unclear what the form of the non- 
Abelian gauge transformations should be on the finite-element lattice. 

Instead, the first foray into non-Abelian finite element gauge theory was based on 
straightforward gauging of the global phase symmetry of the free finite-element Dirac equa- 
tion. The form of the interacting Dirac equation was determined, nonlocally as in the 
Abelian case, and at the same time, the form of the gauge transformations of the vector and 
scalar potentials was determined, in terms of an infinite sequence of nested commutators. 
The Yang- Mills equations were determined analogously. The only thing not explicitly deter- 
mined at the time was the form of the construction of the field strength tensor in terms of the 
potentials; although it was perfectly clear that the process could be continued indefinitely, 
only the first four terms in the sequence in powers of potentials were given. Although 
had been restricted to (1 + 1) dimensions for simplicity, that restriction was easily removed 

i- 

The completion of this construction was only given this spring 0. The essential element 
was the recognition that under the previously-determined gauge transformations, a suitable 
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link operator transformed appropriately. These operators then can be used to transform the 
field strength F^j^ averaged over the finite-element hypercube to the /x, v plane where it may 
be expressed as a path-ordered product of link operators around a plaquette. 

This development makes it possible to revisit Matsuyama's scheme ||^. We will see that 
it is not only possible to formulate a local Dirac equation in the non-Abelian regime, but 
local gauge-covariant Yang-Mills equations as well. The resulting equations are inequivalent 
to those given previously but not so different either, for the previous equations were 
quasilocal, as seen in the simple difference equation given for the interaction terms. But the 
new formulation is still nonlocal in that the field strength that appears in the Yang-Mills 
equation involves the vector potential over the entire previous lattice. 

In the next section we restate the gauge transformation properties of the link operators, 
and give the corresponding construction of the field strength. Then, in Sec. 3 we restate 
Matsuyama's prescription for the Dirac equation, followed by the corresponding local Yang- 
Mills equation. The resulting nonlocality in the construction of the field strength is shown. A 
simple calculation of the axial-vector anomaly in two dimensions is given in Sec. 4. However, 
it is not clear how to extend such calculations to four dimensions because, unlike the original 
formulation in general interactions here break unitarity. In a particular gauge in which 
the transfer matrix is unitary, the current anomalies are computed in the smallest nontrivial 
four-dimensional lattice in Sec. 5. The corresponding calculation in the standard finite- 
element formulation is given in Sec. 6. A discussion of how symmetry breaking occurs here 
is given in the Conclusion and the Appendix. 

II. GAUGE-COVARIANT LINK OPERATORS AND CONSTRUCTION OF THE 

FIELD STRENGTH 

It is convenient to define the link operator using coordinates referring to a given finite 
element: 
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where h is the lattice constant, and the vector potential is an appropriately averaged one: 
The notation here is that m_|_ refers to space-time lattice coordinates other that the one 



/1M 



singled out, m^. Each index in ( |2.1| ) takes on two values, or 1. The result of a detailed, 



constructive calculation is the following simple transformation law for the link operator 

1000 — '^5' ['^'^0000 (-^ A*) 1000 ~ (-^Ai) 1000*^^^1000] ) (2.3) 

where we have assumed that the first coordinate index refers to the /i direction. Then, 
it is easy to see that the "transversely-local" field strength JF^^^ is given by the following 
path-ordered product of link operators around the fi, v "plaquette" : 

where the first index is the /i coordinate and the second index the u coordinate. To construct 
the full field strength, which is forward-averaged over a hypercube with lower left-hand corner 
at m we use the gauge covariant average operators constructed from the links: 

{D.^noooo = \ [{LM^ni{L{)i + {^^o] , (2.5) 

where on the right side we have only displayed the A coordinate (the rest are 0). That is, 
with an overbar representing forward averaging, 

~ '^{.•^m ~l~ ^^ni+l); (2-6) 



the field strength is given byQ 



^The averaging over the finite element on the left side of ( |2.7| ), required by the finite-element 
prescription, is necessary to ensure unitarity. If {F^y)Ya were replaced by (-F^jy)mi even the free 
theory would not be unitary. 



where symmetrical averaging is to be understood. It is immediately obvious that F^^, given 
by ( p.7|) transforms covariantly, 

{5F^^'')^ = ig[5uj^,{F^-)^]. (2.8) 



III. LOCAL FORMULATION OF YANG-MILLS EQUATIONS 

Matsuyama had proposed a local finite-element formulation of a fermion interacting 
with an Abelian gauge field 0. He began by adopting a local form for the fermionic gauge 
transformation, 

(^^m = igSu^ipra, (3.1) 

defined on fields at the lattice sites, rather than in the middle of the finite element as in 
0,^J^. Then covariant derivative and averaging operators can be defined in terms of the 
link operators defined in (|27 



(^A.^)oooo = [(^m)i^i - M ' (3-2a) 

(D^^PUoo = \ [{LMi + M ■ (3-2b) 

Up to an ordering ambiguity then, the gauge covariant Dirac equation is 

J] D^D^ij + J] = 0. (3.3) 



By virtue of ( |2.3| ) this equation ( |3.3| ) is covariant not only under Abelian gauge transforma- 
tions, but under non- Abelian ones as well. Finally, we can transform to a gauge-covariant 
average Dirac field by averaging operators, 

\ A / m 

[t^ 'i/'m using the notation of ( ^.61 )] which transforms covariantly in the sense of 

(5^M = ig^^r^^m- (3.5) 



However, in ||^ it the "free" average of the Dirac field, ipm that transforms covariantly. 

We proceed similarly in the gauge sector. If the field strength like the Dirac field trans- 
formed locally at the lattice sites, 

{6Fnm = ^9[SuJr., {Fnm] (3.6) 

rather than as in gauge covariant derivative and averaging operators could be defined 
by ( p.5|) and 

{D.Fnoooo = I [{LMFni{L{)i - [Fno] , (3.7) 
which yields the local covariant Yang-Mills equation: 

llD^D.F^^'' =f, (3.8) 

where we can adopt the following as a gauge-covariant current (see Sec. 6): 

if)^ = g^^T^^^. (3.9) 

However, ( |3.6|) does not hold! As shown in Sec. 2 the field strength constructed locally in 
terms of the link operators transforms according to ( p.8|) , defined at the center of the finite 
element. However, given F^^, constructed in Sec. 2, we can construct a locally covariant field 
strength /^,^ according to 

iF,.)m= [HDxUA , (3.10) 

which does transform according to ( |3.6| ), and which then satisfies ( p.8|) . Note that, in part, 
( p.lOp undoes the transformation ( p.7|) , so that if we choose an appropriate ordering 

However, as a result of inverting the averaging operators D^, that is solving the difference 
equation ( |3.11D , ff^i, is not local, but depends on vector potentials over the entire prior 



lattice. It appears impossible to have a completely local formalism. 



IV. AXIAL- VECTOR CURRENT ANOMALY IN (1+1) DIMENSIONS 



Let us illustrate the calculational aspects of this scheme in the simplest context, that of 
an Abelian theory in which the only nontrivial element is the Dirac equation (|3.3|) . (We will 
henceforward replace g hy e.) Although the result was stated in 0, it is useful to first revisit 
the two-dimensional case of the Schwinger model, with the fermion mass /i = 0. In terms of 
chiral components, that is, eigenvectors of with eigenvalue equal to ±1, the solution of 
that equation is (for a particular ordering) 

^i:l,n+i = e^^^('^»)-"+^e^^^(-^^)™+^-+^<) , (4.1a) 
= e^^^(-4o)-"+^e-^^('^^)-+^-^tJi,„. (4.1b) 

The current is given by ( p.9|) with T = 1; in terms of chiral components what we wish to 
compute are the following combinations: 

\ '-^t^J I ^ Wm+l,n+l Jm,n Jm,n+1 Jm+l,n/' y±.^cx) 
{ ^fJ-jb ) ~ J^{jm+l,n+l ~ jm,n ~ Jm,n+1 ~^ jrn+l,n) ^ (4.2b) 

where the quotation marks signify finite-element lattice derivatives. We use the solution 
( [4.1| ) to refer all Dirac fields to the intermediate time n + 1 and we evaluate the fermion 
matrix elements according to the Fock space rule |^ 

(V^SVi2g,n) = igy— ^TTT, for M even, (4.3a) 
(^Sr^S,,n) = igy^^^^^, for M odd, (4.3b) 

where g = ±1. In both cases, the vacuum expectation value is taken to zero if m = m'. (As 
noted in |^ the actual value in that case is irrelevant.) Then using ( |4.1| ) in, for example. 



_|_ g-jeh(^o)m,n+i^(+) _|_ g-«e/i(.Ai g-jeh(.Ao)m+l,n+l^(+)^ ^\ (4 4) 
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we easily find that the vector current is conserved exactly^, 

('V) = (4.5) 

while the axial-vector current is anomalous, 

e 1 

= TrrjU- 7T7[sine/i(A)m+2,n+i + sine/i(A)m+i,n+i 

+ sin e/;,[(^o)m+l,n+l — {^o)m+2,n+l ~ (-41)01+2,71] 
— sin e/l[(^o)m+l,n+l + {•^l)m+l,n — ("4o)m,n+l] 

^ . ' , E, {h^O), (4.6) 
M sm vr/M 

where E = "c?o^i — diAo" is the lattice electric field, the familiar finite-element lattice result 

i- 

The above calculation seems quite similar to that given for the nonlocal finite-element 
formulation in and is certainly no simpler. In fact, the calculations are identical if, as 
in P], we choose the gauge Aq = 0. For then the massless Dirac equation can be written in 
terms of the transfer matrix, defined by 

tpn+l = Tipn, (4.7) 

which is to be understood as a matrix equation in the spatial coordinates. Here the transfer 
matrix is 

T = 1±^, (4.8) 

where the covariant derivative, V = —{h/2){Di)~^Di, is defined in terms of (the time 
coordinate n is suppressed) 

(I^i)™,m' = \ (-5™,™' + <5^+i,™'e-^^''(^^)-+i) . (4.9b) 



^In @ results (U) and (0) were established only to 0{h). 
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Equation ( |4.9a|) is inverted as in (17) of P] with the result that the covariant derivative 
operator coincides exactly with that given in (10) of that reference, and hence the same 
conclusions follow. 



V. CURRENT ANOMALIES IN (3+1) DIMENSIONS 



We wish to repeat the above calculation in four dimensions. Again we will set the mass 
H = 0, so we have to solve the following symbolic Dirac equation involving the link operators 
( p.l|) for the simplest possible ordering: 



7°(L3 + 1)(L2 + l)(Li + l)(Lo - 1)^ + ^\Ls + 1)(L2 + l)(Li - l)(Lo + 1)V 
+ 72(L3 + 1)(L2 - l)(Li + l)(Lo + 1)^ + 7'(^3 - 1){L2 + l)iL, + l)(Lo + 1)^' = 0, (5.1) 



where, for example, ij = V'oooo, Lq^j = (Lo)oooi^oooi, and L2Loi) = (L2)oioo(-^o)oioiV^oioi- We 
must solve this system of equations at each lattice site at a given time. For simplicity, let 
us consider the simplest nontrivial rectangular spatial lattice, with the number of sites in 
the 1 direction being Mi = 2, while in the other two directions there is but a single site, 
M2 = M3 = 1. We anticipate the periodic/antiperiodic boundary condition, 



so we expect that the fields should be antiperiodic in the 1 direction. The system of Dirac 
equations reduces to the following simple matrix problem 



V'm+M = (-1) 



(5.2) 



(5.3) 



where the Dirac fields at time 1 and time are 



^1 





(5.4) 



and the matrices are given by 
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and 



^ —Iqc 21qIi —lohc 

Iqc* 2lo lohc* 

—21qIi lolid — /qc? 

\-lohd* kd* 2lo 



\ 



2 c /ic 

-c* -he* 2h 

-/ic/ 2 d 

yhd* -2h -d* 



where 



2/o Iqc /o^ic 
-/oc* -lohc* 2loh 
—lohd 2/o /o*^ 



c- 



7^'5 



-c 2/i -/ic 

c* 2 he* 

-2/i /irf -d 

y-hd* d* 2 J 



(5.5) 



(5.6) 



^0 — (-^0)01, ^0 — (-^0)11, h — {Li)()Q, h — {Li) 



10, 



and 



c = tan y (^3)00 + i tan y (^2)00, d = tan y (^3)10 + z tan y (^2)10, 



(5.7) 



(5.8) 



and we have denoted eigenvalues of ^75 by the ± superscripts, have used the following 
representation of Dirac matrices: 



7°7 = ^750-, 



(5.9) 



and have chosen cxi to be the Pauli matrix. The transfer matrix T is then given by 



T = R-^S. 



(5.10) 



Unfortunately, it turns out in general the theory is not unitary, that is TT"!" 7^ 1. This is 
true even in the temporal gauge, Aq = 0. This is traced back to the failure of the covariant 
derivative, given by (|3.3|) , 



Ao 



(5.11) 
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where symbolically 

Ao = (L3 + 1)(L2 + + 1), Ai = (L3 + 1)(L2 + - 1), 

A2 = (L3 + l)(L2-l)(Li + l), A3 = (L3-l)(L2 + l)(Li + l), (5.12) 

to be skew Hermitian.0 In fact, it is easily verified that although Vi is skew Hermitian, D2 is 
not unless /i/i = 1. So, to proceed, we will simply choose the gauge with Aq = Ai = 0, which 
should suffice for the anomaly calculation. This simplifies the transfer matrix dramatically 
by replacing 



1. 



(5.13) 



Then it is easily seen that T is unitary. (If ip were periodic rather than antiperiodic in the 
1 direction, which would be accomplished by changing the sign of /i in ( |5.5| ) and ( |5.6| ), T 
would not be unitary.) In fact, it will suffice in the following to expand T to bilinears in c 
and d. It then turns out to be 

( „, A 1 I t n \ 



V 


d 







-c* 


V 







1 + e 





V* 


c 





-1 + ^ 


-d* 


V* 



\ 



where we have abbreviated 



1 



^ = -dd*, = — cc*, V = — dc*. 



V — 

d* V 

-1+^ V 

l + i c* 

1 



c 1 + e 
-1 + ^ 
-d 



V 



(5.14) 



(5.15) 



In the following we will use the block form of T, which refers to specific spatial coordinates: 

T= - - ^ (5jg) 



^10 2^11 



^It is easily seen that no other ordering will resolve this problem. For the given ordering, 1)2 = 
(Li + 1)~^(L2 — 1)^^{L2 + ^){Li — 1); the inner factor involving L2 is clearly skew Hermitian, but 
the appearance of the Li terms destroys that property. 
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We now take matrix elements of the current in Fock-space states defined in terms of the 
momentum-space expansion of the free Dirac field 

^m^i = T.S Ku.se'^'^'/'^"''''/'' + 4.t;,.e-'(^+^/^)-^-/^0 • (5-17) 

s,p 

Because of the particular choice of gauge (|5.13|) the current is proportional to the usual 
finite-element one, and the matrix element of for example, is 

ifrnj = C^{[{4il + T))i + (^1(1 + T))o][((l + Tt)^Oi + ((1 + Tt)^i)o]), (5.18) 

where the Dirac fields are those given in ( ^.41) , 

C = cos2y(A)ooCos2^(^3)oo, (5.19) 

and the matrix subscripts refer to the spatial coordinate. The matrix elements are evaluated 
using the following easily derived formula for M = 2: 

ii^lTijm') = ^{Smm'tlT + 2l 6 mm' tl -f'^^'V) . (5.20) 

Then a straightforward calculation gives the result 

(An) = ;^(l + A)C, (5.21) 

where 

A = i {dd* - cd* - dc* -cc*). (5.22) 

A similar calculation reveals that (j^) = 0. The axial- vector current matrix elements vanish. 
Consequently, for this tiny lattice, there is no axial-vector current anomaly, but there is a 
vector current anomaly, 

= -^{[(^2)?i + (A)ll(A)01 + (A)^i - (A)?o - (A)l0(^2)00 - (^2)^0] 

+ [2^3]}, (5.23) 
a somewhat curious result which will be discussed below. 
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VI. FINITE-ELEMENT CURRENT ANOMALY CALCULATION 



Because of the various difficulties seen in the calculation exhibited in the previous sec- 
tion, we return now to original formulation. There, the unitary transfer matrix in (3+1) 
dimensions is, in the gauge = 0, |0] 

1 + 7°7 • D 



T 

1-707-1)' 
where (the time coordinate n is suppressed) 



(6.1) 



'( 1) ' ' (^mim' COS (^miTTt' ^ sin (^^^^' ) sec (^^ 



m I m I ) 



(6.2) 



where (and now all the local dependence on the other spatial coordinates is also suppressed) 

M 



eh 



(6.3) 



mi=l 



and 



M 



Cm„m^ = X] {m'l - rui) sgn {m'l - m•)C„^'/ 
Here, the sign function is 



(6.4) 



m '=1 



sgn(mj 



1, m > 0, 



-1, m < 0. 



(6.5) 



Again, let us consider the smallest possible nontrivial lattice, with the number of lattice 
points in the 1 direction being Mi = 2, while the 2 and 3 directions have but one site, 
M2 = M3 = 1. This leads again to an 4 x 4 transfer matrix for each chirality, which for 



small eh A is, in the block form given in ( 5.16 ) 



± 
00 



A 



5 -^01 ~ -^10 



a P 
f3* -a 



\ 
I 



T+ — tA — _ 

^10 ~ ^01 ~ 2 



a p 
(3* -a* 



^ 11 



-A ±^ 
+^ -A 



(6.6a) 



(6.6b) 
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Here A is defined above in ( |5.22 



~X = ^(^dd* + cd* + dc* -cc*), (6.7) 



c and d are as in ( p78 



1 

a = -(4 + 42ti - Ati - dd* + d*c - c*d - cc*), (6.8a) 

a = i(-4 - 4iti + 4ti + rfrf* + d*c - c*d + cc*), (6.8b) 

f3 = = d{l + iti) - c{l -iti), (6.8c) 

/3 = = d(l-iti) -c(l + iti), (6.8d) 



and 



ti = tanC^^^ s = -e^^'^i-^^) secC^^^ (6.9) 

Even though the transfer matrices are quite different (the form in ( |6.6|) is more comphcated 
because Ai 0), the results of the calculation are very similar: The axial- vector anomaly 
is zero, while the vector anomaly is 

(' V) = -^[(^2)21(^2)11 + (^3)21(^3)11 - (^2)20(^2)10 - (^3)20(^3)10]. (6.10) 



Here, the current is not ( ^.9] ) but 

Jm,n = e^m,n7^^m,n, (6.11) 



the averaging being over the finite element without the link operators. The result ( |6.10[ ) 



which is manifestly gauge invariant, appears to be, like ( |5.23D , a lattice artifact.^ It is not a 
lattice version of F"^, as one might anticipate. However, we should not be discouraged, since 
the calculation given is for a truly tiny, unrealistic lattice. The extension of this calculation 
to larger lattices will be presented elsewhere. 



^It is probable that this result reflects the rectangular nature of the lattice considered here. A 
similar "anomalous" anomaly was found when the space-time lattice was not chosen to be square: 
That is, when the lattice spacing in the time direction was not equal to that in the space directions. 
See @. 
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VII. HERMITICITY, UNITARITY, AND THE EXISTENCE OF A SYMMETRY 

CURRENT 



In this paper we pursued a variation on the finite-element formulation of lattice gauge 
theory which seemed at first sight very promising. The idea was to use the link variables, 
in term of which the local field strength was constructed, to express the Dirac and Yang- 
Mills equations in local form, rather than the form involving the entire prior lattice given 
previously. A locally gauge-invariant construction can indeed be done, one which is inequiv- 
alent to the earlier formulation. However, the new scheme is less advantageous than it first 
seemed, and ultimately fails to be consistent: 

• The field strength which appears in the Yang-Mills equation is not locally constructed 
in terms of local link operators. It appears impossible to have a local formulation 
consistent with local non-Abelian gauge symmetry. 

• Moreover, detailed calculations, even in the Abelian theory where the equations are 
local, turn out to be no simpler, and perhaps more complicated than those in the 
original formulation. 

• Disastrously, unitarity is violated. Explicitly, we have seen in the Abelian case that the 
transfer matrix is not unitary, even in a temporal gauge. This is traced to the fact that 
the covariant derivative operators are not skew Hermitian. In contrast, the original 
formulation is manifestly unitary (canonical). It should be recalled that preservation 
of the canonical commutation relations at the lattice sites was the original motivation 
for adopting the finite-element prescription for field theory on a Minkowski lattice. 

We can interpret these results in a positive light. It was always apparent that, although 
some arbitrary choices had to be made to implement local gauge invariance, the require- 
ment of that invariance, that is, that the transformation equations could be "integrated," 
was rather rigid. The findings presented here strengthen that conclusion. The fact that 
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we have now shown that a rather natural and attractive alternative gauging process is ulti- 
mately inconsistent makes the pursuit of extracting physical information from the orginal, 
consistent, approach that much more compelling. 

Finally, we should add some remarks about the current employed in both formulations 
studied here. It should be noted that the choice of current was essentially arbitrary, subject 
only to the requirement that it be locally gauge covariant. It is essential to note that the 
current cannot be derived from the Dirac equation. In fact, because our equations of motion 
cannot be derived from an action, there is no connection between symmetry (say chiral 
symmetry) and conservation laws (say axial- vector current conservation). In our consistent 
formulation this is because the Dirac equation is asymmetric between past and future. If one 
were interested in a Euclidean formulation, one would choose a Dirac equation symmetric 
in the fourth coordinate, and an action, and corresponding current could be constructed. 
That current would, however, also involve all lattice sites in the fourth coordinate, and 
therefore would be be unusable in the Minkowski context, where we wish to solve the operator 
equations of motion by time-stepping through the lattice. This Euclidean construction and 
its Minkowski failure is sketched in the Appendix. 
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APPENDIX: CURRENT CONSTRUCTED FROM EUCLIDEAN LAGRANGIAN 



In the text we simply assumed a form of the current ( ^.9] ) and ( |6.11| ) which was manifestly 
gauge invariant. We are at liberty to do so, because the Minkowski finite-element equations 
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of motion are not derivable from a Lagrangian. The current cannot be derived from the Dirac 
equation, but is an independent source for the Maxwell equations, see Ref. 0,^. However, 
if one were to work in Euclidean space-time (periodic or antiperiodic in all four directions), 
it is possible to construct an action from which the equations of motion are derivable, and 
which therefore supplies a lattice current. The fermion part of that action is (a factor of i 
is absorbed in going to Euclidean space) 

w^/ = ^'E^^(rr-^ + '^/^) r = (7V,7°)- (ai) 

m,m' ^ ^ m.m' 

If we vary ( |Al| ) with respect to ip'^ we obtain the Euclidean Dirac equation 

■V + if?j^ = 0. (A2) 

where, as in (|Al|), a four-dimensional scalar product is implied. 

Given an action, we can construct a conserved vector current by making a local gauge 
transformation 

Sipm = ieSQjrii'm- (A3) 

Because the Dirac equation, and hence the action, is invariant under the global version of 
(|X3|), 6Qrn = 6Q = constant, we must have by the action principle 

SWf = = -h^ Yl Jinl{S^m,,m^ - 5^^„^.-l,mJ, (A4) 
m 

from which we read off the conserved current 

= -e ^ i^ln'^,mj'^m'^,m^ sgn(mi - m[) sgn(mi - 

m'- ,m'' 

x(-l)-H<secC(^)exp(-^e„,,^.a,,^»)• (A5) 

(The same result, of course, can be obtained by varying V ( |6.2| ) with respect to .4.^. j„^.) 
The expression for this Euclidean current has been simplified by deleting constant terms. It 
is easy to verify explicitly that this current is both conserved and gauge invariant. Similarly, 
by making a chiral transformation, 
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Sipm = 7°75'5^mV'm, (A6) 



we can construct the axial- vector current J^^, which has the form of (|A5|) with the replace- 
ment 

er^j^j,r^Tl Tl = i-t^,Y,-i75). (A7) 

By construction, these currents possess no anomalies. However, they appear to be completely 
unacceptable, because they are horribly nonlocal. In particular, they possess no Minkowski 
analogues, in the sense that it is not possible to analytically continue back to real unbounded 
times. Crucial to our formulation is the propagation of the operators from past times to 
the present time, so that we can solve for the field operators by time-stepping through the 
lattice. The Euclidean current ([A5|) involves fermion field operators at all Euclidean times, 
which would make it impossible to solve for the operators at time n in terms of operators at 
earlier times. Therefore, for the considerations of the text we use the gauge-invariant current 
( p.9|) and ( |6.11|) and their axial analogues, currents which can and do possess anomalies. 



It is further illuminating to note that if we were to use the current ( |A5[ ) in a one-loop 
lattice calculation of the vacuum polarization in two dimensions, we would find a vanishing 
anomaly, rather than the value e^/vr reported in This is because species doublers occur 
in the action defined by (|A1|). Such doublers are absent in the finite-element scheme based 
on equations of motion and the current ( |6.11D . 
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